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Quantum information requires protection from the adverse affects of decoherence and noise. This review 
provides an introduction to the theory of decoherence-free subspaces, noiseless subsystems, and dynamical 
decoupling. It addresses quantum information preservation as well as protected computation. 



I. INTRODUCTION 



The protection of quantum information is a central task in quantum information processing [ 1 1. Decoherence and noise are 
obstacles which must be overcome and managed before large scale quantum computers can be built. This review provides an 
introduction to the theory of decoherence-free subspaces, noiseless subsystems, and dynamical decoupling, among the key tools 
in the arsenal of decoherence mitigation strategies. It is based on lectures given by the author at the University of Southern 
California as part of a graduate course on quantum error correction, and as such is not meant to be a comprehensive review, nor 
to supply an exhaustive list of references. Rather, the goal is to get the reader quickly up to speed on a subset of key topics in the 
field of quantum noise avoidance and suppression. For previous reviews overlapping with some of the theoretical topics covered 
here see, e.g., Refs. J2][3l. 

The review is structured as follows. Section|lI]introduces decoherence-free subspaces (DFSs). Section|lII]defines and analyzes 
the collective dephasing model, and explains how to combine the corresponding DFS encoding with universal quantum compu- 
tation. Section IV considers the same problem in the context of the more general collective decoherence model, where general 
noise afflicts all qubits simultaneously. SectionJV] introduces and analyzes noiseless subsystems (NSs), a key generalization of 
DFSs which underlies all known methods of quantum information protection. The NS structure is illustrated with the three- 
qubit code against collective decoherence, including computation over this code. We then proceed to dynamical decoupling 
(DD). Section VI introduces the topic by analyzing the protection of a single qubit against pure dephasing and against general 
decoherence, using both ideal (zero-width) and real (finite-width) pulses. Section [VIJ| briefly discusses DD as a symmetrization 
procedure. Section [VIII| discusses combining DD with DFS in the case of two qubits. Section [IX| addresses concatenated dy- 
namical decoupling (CDD), a method to achieve high-order decoupling. In the final technical Section |Xj we come full circle 
and connect dynamical decoupling to the representation theory ideas underlying noiseless subsystems theory, thus presenting a 



unified view of the the approaches. Concluding remarks and additional literature entries are presented in Section XI 



II. DECOHERENCE-FREE SUBSPACES 



Let us begin by assuming we have two systems: A and B, defined by the Hilbert spaces Ha and Hb, respectively. In general, 
the dynamics of these two systems are governed by 

H = H A + H B + H AB , (1) 

where Ha and Hb correspond to the pure dynamics of system A and B, respectively, and Hab is the interaction between the 
two systems. For simplicity assume Ha = and Hb = 0, and let us focus on the dynamics of A resulting from the interaction 
Hamiltonian. Recall that using the Kraus representation Q), we can effectively study the reduced dynamics of A for an initial 
state pa(0), where 

PA {V)" PA(t) = Y J K a {t) PA {0)Kl(t) (2) 

a 

after the partial trace over system B is completed. The Kraus operators K a (t) satisfy the relation Y, a Ka(t)Ka(t) = I A Vt, 
where I a is the Identity operator on system A. In general, we see that the Kraus operators result in non-unitary evolution in the 
system A Hilbert space. Therefore, let us define decoherence as follows: 

Definition 1. A open system undergoes decoherence if its evolution is not unitary. Conversely, an open system which undergoes 
purely unitary evolution is said to be decoherence-free. 

It is from this definition of decoherence that a DFS will be defined for system A. 

Since the DFS defines a portion of the system Hilbert space where decoherence is not present, let us partition the Hilbert 
space into two subspaces Ha = Hg ®Hn- The section of the Hilbert space not effected is denoted by Hg, the "good" portion, 
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and TL n denotes the decoherence-effected Hilbert space. From the definition of decoherence we find that if Ha is to remain 
uneffected it must evolve unitarily. Therefore, assume that it is possible to partition the Kraus operators as 

K a {t)=g a U®B a (3) 

such that U defines a unitary operator acting on Hg, 9a 6 C, and B a is an arbitrary (possibly non-unitary) operator acting solely 
on "Hat. Note that by the completeness relation of the Kraus operators, 

£| 5q | 2 = 1, Y, B 1 B » = In- (4) 

a a 

Assume that the initial state is partitioned as well, as 

Pa(0) = Pg (0)®Pn(0). (5) 
In terms of the newly partitioned Kraus operators Eq. Q then becomes 

P A -P' A -T(9 a ue B a ) PA (0) ( 9 * a lP e Bl ) = ( U <*W* Zq B J N{Q)Bi ) . (6) 

The remarkable thing to notice about this last result is that pc evolves purely unitarily, i.e., it satisfies the definition of 
decoherence-freeness. Thus, we conclude that if the Kraus operators have a block diagonal form as in Eq. Q then a DFS 
Hg exists. 



A. Deutsch's Algorithm 

As a first example of error avoidance utilizing the DFS construction, we can consider the first known algorithm that offers a 
quantum speed-up. Deutsch's Algorithm [ 1 1 presents a simple decision problem in which the goal is to decide whether a function 
is constant or balanced. Let f(x) : {0, l} n •-»■ {0, 1} (in decimal notation x e {0, 1, . . . , 2™ - 1}), denote the function, where if 



/(*) 



0, Vxor 

1, Vi 



(7) 



then the function is called constant and if 



fix) 



0, half the inputs 

1 , half the inputs 



(8) 



then the function is called balanced. Classically, we find that in the worst case, making a decision on whether f(x) is constant 
or balanced requires a minimum of 2 n /2 + 1 total queries to /. Deutsch and Jozsa found that the exponential cost in /-queries is 
drastically reduced by considering a quantum version of the algorithm. 

The equivalent version of the decision problem can be recast in terms of the quantum circuit: 



11} 



H 



1^2} 



1^3} 



1^4} 



where each classical bit n corresponds to a qubit. The unitary operator U performs the query on f{x) by 

Uf ■ \x) \y) h> \x) \y e f(x)) (addition mod 2) (9) 

and H represents the Hadamard gate: H\0) = \+) and H \l) = |-}, where |±) = (|0) ± |l})/\/2. 

In order to illustrate Deutsch's algorithm for the quantum circuit above, consider the single qubit version (n = 1). In this 
case there are four functions, two of which are constant and two of which are balanced: {/o(a;) = 0, fi(x) = 1} (constant), 
if2(x) = x,fz{x) = x} (balanced), where the bar denotes bit negation. Clearly two classical queries to / are required to tell 
whether / is constant or balanced. 



3 



Initially the total system state is given by \ipi) = ]0) |1). Applying the Hadamard gate, the system state becomes |-02 ) = | + ) | - )- 
Applying the unitary operator Ut, the resulting state is 



\*h) = u f \+) \-) = u f 
1 



(|00>-|01> + |10>-|11» 



2 (|0, e /(0)) - |o, l e /(0)) + |1, e /(l)) -|l,le f(i))) 
i(|0,/(0))-|0,/(0)) + |l,/(l))-|l, /(!))). 



(10a) 
(10b) 
(10c) 



Applying each constant and balanced function to IV's)* we find 



|V>3> 



{ fo- 


|00) 


-|01) + 


10}- 


111)=+ 


+-> 


fl- 


|01) 


- |00) + 


11}- 


|10) = - 


+-> 


/ 2 = 


|00) 


-|oi) + 


11)- 


|10) = + 


--) 


/3 = 


|01) 


-|00) + 


10)- 


111) = - 


--) 



(11) 



for fj(x) € {0, l,x, x} as defined by Eqs. ^ and (|8j. The remaining Hadamard gate yields the final state 



r/o 

h 
h 

1/3 



+ |0-) 

-|o-) 

+ ll-> 
-IH 



(12) 



and the characteristic of the function is determined by measuring the 1st qubit: a result of indicates a constant function, a result 
of 1 a balanced function. Thus, remarkably, we find that the quantum version only requires a single query to the function, while 
the classical case requires two queries (this scenario is the original Deutsch algorithm). 

The circuit depicted above can be subjected to a similar analysis in the n-qubit case (the Deutsch-Jozsa algorithm) and the 
conclusion is that the quantum version of the algorithm still requires only a single /-query, thus resulting in an exponential 
speed-up relative to its classical counterpart in the worst case. 

B. Deutsch's Algorithm With Decoherence 

To gain an understanding of how a DFS works we can look at the Deutsch Problem with added decoherence. We can consider 
the circuit diagram for the single qubit Deutsch algorithm, but introduce a dephasing element as follows, where the dotted box 
denotes dephasing on the top qubit only: 



10) 
|1> 



H 



H 



Z 



HZ} 



92 



P2 



P3 



H 



Pa 



The Kraus Operators governing the dephasing of p 2 are: 

K a = y/l-ph ® h 
K x = X /]5Z 1 ® I 2 



(13a) 
(13b) 



With probability (1 - p) nothing happens. However with probability p the first qubit experiences dephasing. If we follow 
the density matrix states through the algorithm we can see the effect this dephasing has on our result. As before we have 
px = |01) (01| and, p 2 = |h — ) (h — |. By applying the Kraus operators we find the state after dephasing to be 



p' 2 = K p 2 Kl + Kip 2 K{ 

= (i-p)p 2 + P (z®i)\+-){+-\(z®iy 

It is easy to check that Z\ + ) = |-), thus we find 

p2 = (l-p)\+~){+~\+P\—){—\- 



(14a) 
(14b) 

(15) 
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Using Eq. |9} we can compute ,03 



P3 



(1 - p) |+-> <+-| + p I > < — 

(l-p)|_)(_|+p|+-><+- 



/o,/i (constant) 
f 2 ,h (balanced) 



After applying the final Hadamard we see that the state we will measure is 



Pi 



(1 - p) |0-) (0-1 + p |1-) (1-1 f , h (constant) 
(1 -p) |1-) (1-1 +p|0-) (0-1 f 2 J 3 (balanced) 



(16) 



(17) 



If we now measure the first qubit to determine whether the function is constant or balanced, with probability p we will misidentify 
the outcome. For example, if we obtain the outcome 1, with probability p this could have come from the constant case. But, 
according to Eq. ([12} the outcome 1 belongs to the balanced case. 

It is possible to overcome this problem by use of a DFS. Let us again modify the original circuit design. We can add a third 
qubit and then define logical bits and gates. 



\0l) -F- 

11} 



Pi 















ZZ - 







H 



u 



fL 



P2 



P-2 



P3 



Pi 



Now the Z dephasing acts simultaneously on both top qubits, that comprise the logical qubit \0l)- In this case the Kraus 
operators are 



K = y/l^pIII, Ki = ^ZZI , 
where ZZI = a z <g> a z <g> /. Recall the requirements for a DFS. The Kraus operators, as in Eq. ([3]), must be of the form 



( g a U \ \ 

{ \B a ) 



(18) 



(19) 



and the state must be initialized in a good subspace, i.e., ps = pa © Pn> where the direct sum reflects the same block structure 
as in Eq. ( |19) , If these conditions are met then p' s = Y, K a psK^ a = UpaU^ © p' N . In other words the evolution of pq is 
entirely unitary. We start by checking the matrix form of the Kraus operators. K a is simply the identity and trivially satisfies 
this condition. We can check the ZZ portion of K\ since that is what will act on our logical qubit. 



ZZ : 



(1 0\ 
0-100 

0-10 

\0 1/ 



00 
01 
10 

11 



(20) 



This obviously doesn't fit the required matrix format, in that there is no block of l's like in Kq. However with a simple reordering 
of the basis states we obtain the following matrix 



ZZ : 



/l \ 

10 

0-10 

V0 -1/ 



00 

11 

01 
10 



(21) 



Now we have a 2 x 2 matrix of l's, so both ZZ and the identity matrix act as the same unitary on the subspace spanned by |00) 
and 1 1 1 ) . The full matrix, K±, then takes the form 



Ki = VP 



fh 








0\ 


OOz 





h 








lli 








-h 





Oli 


V0 








hi 


lOi 



(22) 
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where I 2 denotes the 2x2 identity matrix, i = {0, 1}, and Kq is the 8x8 identity matrix. Thus we see that both Kq and K\ have 
the same upper block format, namely U = I4 X 4, where g = \A ~P an d 9i = \fV- Now we can define our logical bits |0l) = |00) 
and \\ L ) = 1 1 1 ) . With these states we can construct our logical Hadamard. 



(Hl) 



4x4 ' 



H 








V 



(23) 



Here the logical Hadamard acts as a regular Hadamard on our logical qubits. 

H L \0 L ) = \+ L ) (24a) 
H l \Il) = \-l) (24b) 

where \±l) = ± The other unitary action of Hl, namely V, we don't care about. Similarly we can construct a 

logical Uf. 



\ w ) 



(UfL) 8x8 = \-^—\. (25) 



Again, Ujl acts as Uf on our logical bits, and W we don't care about. Neither V nor W affect our logical qubits in any way. 
Now that we have set our system up we can apply the Deutsch algorithm again to see if the DFS corrects the possibility of 
misidentifying the result. Our system begins in the state pi = |0l)|1)(0l|(1| and after applying the logical Hadamard we get 
P2 = \+l) \~) (~\- Now we can apply the Kraus operators to see the effect of the decoherence. K has no effect other than 
to multiply the state by \J\--p because it is proportional to the identity matrix. It is enough to examine the effect of K\ on the 
state \+l)- 

K 1 \+l) = -j=(K 1 \0l) + K 1 \I l )) (26a) 

= -j=y/p(I\0 L ) + I\l L )) (26b) 

= Vp\+l) (26c) 

Therefore p' 2 = p 2 . The decoherence has no effect on our system and the rest of the algorithm will proceed without any possibility 
of error in the end. 



C. A Classical Example 

We can look at a classical case to better understand this effect. Assume we have three parties: Alice, Bob, and Eve. Alice 
wants to send a message to Bob, and Eve wants to mess that message up. Let's also assume that the only way in which Eve can 
act to mess up the message is by, with some probability, flipping all of the bits of the message. If Alice were to send only one 
bit to Bob there would be no way of knowing if that bit had been flipped. But let's say Alice is smarter than Eve and decides to 
send two bits. She also communicates with Bob before hand and tells him that if he receives a 00 or 1 1 he should treat it as a 
logical 0, and if he receives a 01 or 10 to treat it as a logical 1. If this scheme is used, Eve's ability to flip both bits has no effect 
on their ability to communicate. 

This is an example of parity conservation. The logical is even parity and the logical 1 is odd parity. Encoding logical bits in 
parity in this way effectively hides the information from Eve's bit flip error. 



D. Hamiltonian DFS 



Discussing DFS in terms of Kraus operators works well, but we'd like to develop a bottom-up understanding of the DFS 
concept, using Hamiltonian evolution. Assume we are given a system in which our computation is occurring, and a bath that is 
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connected to the system. The Hamiltonian governing the whole system can be written as usual as 

H = H s ®I b + I s ®H b + Hsb, (28) 

where H$ acts only on the system we are interested in, Hb acts only on the bath, and H$ b governs the interaction between the 
two. Assume also, without loss of generality, that the interaction Hamiltonian can be written as 

H SB = Y,S a ®B ai (29) 

a 

where each S a is a pure-system operator and each B a is a pure-bath operator. The Hilbert space can be written % = T~Ls ® T~Lb, 
and Hs = He ® % where 

H G = Span{| 7l )} (30a) 
H N = Span{|i/ fc )} (30b) 
H B =Span{|&)} (30c) 

To formulate a theorem we also need the following assumptions: 

1. The system state is initialized in the good subspace: 

PS = Pg = E r u"l7*)(7jl (31) 

i,j 

2. The basis states of the good subspace are eigenvectors of the interaction Hamiltonian: 

S a |7i> = c a |7i) ,c Q e C (32) 

3. The basis states of the good subspace, when acted on by the system Hamiltonian, remain in the good subspace: 

H s \ 1% )eU G (33) 

With these assumptions in hand we can posit the following theorem. 

Theorem 1. Assuming 1-3, the evolution of the system, ps, can be written as 

Tr B [U ( Ps (0) 9 p B (0)) U f ] = Uspa (0) U f s (34a) 
Ush}en G (34b) 

where U = U (t) = e^'. 

Proof. Using equations ( |29| ) and ( [32] > we can write 

(Is ®H B + H sb ) hi) s ® \Pj) B = hi) ® (Hb Ift )) + E c « l7i) ® B « |ft ) (35a) 

a 

= \li)®\Y J c a B a + HB\\P j ) (35b) 
= | 7i )®F c |^>, (35c) 



where He acts only on the bath. Applying this to equation ( |28) we find that the complete Hamiltonian can be decomposed into 
a portion that acts only on the system and a portion that acts only on the bath. 

H \n) ® ) = (H s ®I B +Is® H c ) |7i) ® \Pj) . (36) 
If we plug this form of the Hamiltonian into the unitary evolution matrix we get 

U (t) = e -^ H s®iB+i 3 <&H c )t = 9 (37) 

where U x (t) = exp(-itH x ), x = S,C. 

To find p(t) we apply this unitary to p(0) with ps(0) = Pg(0). 



p(t) = U(p s (0)®p B (0))U T (38a) 
Jl®Ucp B (Q)U c 



UsPg(0)UI®UcPb(0)U c (38b) 
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We find the state of our system of interest ps(t) by taking the partial trace. 

Ps(t) = Tr B [Usp G (0)Ul U cPB (0)U f c ] (39a) 
= U sPG (0)Ul (39b) 

□ 

Thus Theorem [T] guarantees that if its conditions are satisfied, a state initialized in the DFS will evolve unitarily. 



III. COLLECTIVE DEPHASING 



A. The model 



Consider the example of a spin-boson Hamiltonian. In this example, the system of qubits could be the spins of N electrons 
trapped in the periodic potential well of a crystalline lattice. The bath is the phonons of the crystal (its vibrational modes). We 
also assume that the system-bath interaction has some symmetry in the sense that the interaction between the spins and phonons 
is the same for all spins, e.g., because the phonon wavelength is long compared to the spacing between spins. This assumption 
is crucial for our purpose of demonstrating the appearance of a DFS. If the potential wells are deep enough then the motional 
degrees of freedom of the electrons can be ignored. Let i denote the index for the set of N electrons in the system (the same as 
the index for the set of occupied potential wells in the solid), let k denote the vibrational mode index, b\ \n) = \/n + 1 \n + 1) is 
the action of the creation operator for mode A; on a Fock state with occupation number n, bk \n) = ^/n \n - 1) is the action of the 
annihilation operator for mode k. With of the Pauli-z spin operator acting on the i th spin, the system-bath Hamiltonian is 

HsB = Zgl, k <T*®(b k + bi). (40) 

The symmetry assumption implies 

9i,k=9k, (41) 
i.e., the coupling constants do not depend on the qubit index. The system-bath Hamiltonian can then we written 

H SB = Y,^®Y,9k-(h + bi) (42a) 

i k 

= S Z ®B Z , (42b) 

where 



S z = ^af, B z ^Y,9 z k (bk + bl). (43) 

i k 

If these conditions are met, the bath acts identically on all qubits and system-bath Hamiltonian is invariant under permutations 
of the qubits' order. The operator S z is a collective spin operator. 



B. The DFS 

Now consider the case of N = 2 in light of the DFS condition Eq. ( |32) : 

N = 2^- S z = Z®I + I®Z. (44) 

Thus 



|00> 




2 - 10 


0) 


^c z = 2 


|01) 




|01) 


- |oi> = o 


^c z = 


|00> 


s. 


|io> 


-|10) = 


^-c z = 


|00) 




-2- 


111} 


^c z = -2 



(45) 



8 



It follows that the DFS's for the two spins are: 



H N=2 (2) 
Hat =2 (-2) 



{|00)} 

{|oi),|io>} 

{111)} 



(46) 



where we used the notation T~Ln(c z ) to denote the "good" subspace Hq for N qubits, with eigenvalue c z . 
In the Hn=2(0) DFS, there are two states, so we have an encoded qubit: 



For three spins we have: 



|0) = |01) logical 
|T) = 1 10) logical 1 



N = 3 => S z = Z ®J <g> / + I® Z®I + Z 



(47a) 
(47b) 



(48) 



Similarly, it follows that the DFS's for the three spins are: 

Hiv= 3 (3) = 

Wjv=3(-1) = 

^7V=3(~3) = 



{|000)} 
{|001),|010),|100» 
{|011>, 1101) ,1110)} 

{llll)} 



(49) 



We find that there are two possible encoded qutrits for N = 3, one in Hjv=3(1) and the other in Hat=3(-1)- 

In general, the DFS "Hat(A) is the eigenspace of each eigenvalue (A = c z ) of S z . It is easy to see that the number of spin-ups 

(0's) and the number of spin-downs (l's) in each eigenstate is constant throughout a given eigenspace. This corresponds to the 

value of total spin projection along z. In fact, for arbitrary N, 



c z = #0 - #1 



(50) 



Figure [T] known as the Bratteli Diagram, shows the eigenvalues (y-axis) of S z as N (x-axis) increases. It represents the 
constellation of DF subspaces in the parameter-space, (AT, A). Each intersection point in the figure represents a DFS. Each 
upward stroke on the diagram indicates the addition of one new spin-up particle, |0), to the system, while each downward stroke 
indicates the addition of a new spin-down particle, |1), to the system. From this we see that the dimension of the eigenspace with 
eigenvalue A is given by the relation: 



&m{H N {c z )) = 



(51) 



The highest dimensional DFS for each TV is thus 



max{ dim( H jy ( c 2 ) ) } 



(1) TVeven 
(> +1 ) TVodd 



Given a D-dimensional DFS, H: 



# of DFS qubits in U = log 2 D. 
This formula should be thought of as giving the rate of the DFS code, i.e., 

_ # of DFS qubits in U log 2 D 



# of physical qubits 



N 



(52) 



(53) 



(54) 



For a DFS of given dimension it may be preferable to think in terms of qudits rather than qubits. For example, for the DFS 
%3(-l), the dimension D = 3, and so this DFS encodes one qutrit. 

Any superposition of states in the same DFS/eigenspace will remain in that particular DFS. But a superposition of states in 
different eigenspaces will not evolve in a decoherence-free manner. 
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12 3 4 

N 



FIG. 1. Bratteli diagram showing DFSs for each N. Each intersection point corresponds to a DFS. The number of paths leading to a given 
point is the dimension of the corresponding DFS. The diagram shows, e.g., that there is a 2-dimensional DFS for N = 2 at c z = 0, yielding one 
encoded qubit. 



C. Universal Encoded Quantum Computation 

From here on, encoded qubits will be called 'logical qubits'. To perform arbitrary single qubit operations, we need to be able 
to apply any two of the Pauli operators on the logical qubits. On our example system, 7^2(0), we have: 

|0)e|01) 
|1)e]10> 

\(j>) : DFS encoded logical qubit, a |0) + b |1) 
U : Logical operator on the DFS qubits 



Thus the logical Pauli-z operator is 



Z\0) = |6)1 y_ 7 ^ T ^\ |01) Z - 101) = |0) 



Z =Z®I^{ !" [ z%1 !" ' z. . (55) 



z\l) = -|i) J llio) - -|10)= -|1) 

and the logical Pauli-a; operator is 

•v ") = Dl .. v Vv v ..( "D VvX |io)= i) 



. X = X®X=>{ [ '-x 1 ' ' ' . (56) 
x\i) = |0)J \|10)^ X |01) = |0) 

In general, suppose we have 2N physical qubits all experiencing collective dephasing. We can pair them into N logical qubits, 
each in ^(O), and perform Z or X logical operations on the i th logical qubit using the following operators: 

Z i = Z 2i -i®hi, (57a) 
X i =X 2i - 1 ®X 2i . (57b) 
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We can define arbitrary rotations about the logical X or Z axis as = exp[iX6] and Rz(4>) = exp[iZ(f\, An arbitrary 

single logical qubit rotation (an arbitrary element of SU(2)) can then be obtained using the Euler angle formula, as a product of 
three rotations: Rx(9 2 )Rz{4>)Rx{Qi). 

To generate arbitrary operators on multiple qubits, we need to add another gate to the generating set: the controlled phase 
gate. The logical controlled phase gate can be generated from Z% ® Zj = Zi® Z 2 j-i. Thus a Hamiltonian of the form 

H s = Y^uzMZi + Z^xMXi ■ Y.n,,(l)X, ® Zj (58) 

i i 

not only does not take the encoded information outside the DFS ^2(0) of each of the N encoded qubits, i.e., satisfies the DFS 
preservation condition Eq. (33) , it is also sufficient to generate a universal set of logical gates over the logical DFS qubits. More- 
over, this Hamiltonian is composed entirely of one- and two-body physical qubit operators, so it is physically implementable. 



IV. COLLECTIVE DECOHERENCE AND DECOHERENCE FREE SUBSPACES 



The collective dephasing model can be readily modified to give the more general collective decoherence model. The interac- 
tion Hamiltonian has the following form: 

Hsb = £ [flf, fc of ® (b k + b\) + gl k at ®b k + gl k a; ® b{] (59) 

where 

^^t 17 ^^) (6°) 

corresponds to the raising (+) and lowering (-) operators for the spin of a particle respectively, i.e., 

a + |0) = |l) (61a) 

cr + |l) = (61b) 

a~ |0) = (61c) 

a-|l) = |0) (61d) 

where here corresponds to the null vector and should not be confused with the |0) state. Thus, e.g., a + = |1)(0|. It should 
be noted that the first term in the summation of ( |59) corresponds to an energy conserving (dephasing) term while the second 
and third terms correspond to energy exchange via, respectively, phonon absorption/spin excitation, and spin relaxation/phonon 
emission. 

By assuming that all qubits are coupled to the same bath (thereby invoking permutation symmetry): 

9 r k = 9k, Vfc, ae{+,-,z} (62a) 

B + B_ B z 

^H SB = Y, o\ ® £ 9 + k h + E < ® E 9k b\ + E ff '®E5fe (bk + b\ ) (62b) 
= E S a ®B a , (62c) 

ae{ + ,-,z} 

where 

Sa = f>? ( 63 ) 

i=l 

is the total spin operator acting on the entire system of N physical qubits. We can derive the following relations directly from 
the commutation relations of the Pauli matrices: 



[S ± ,S Z ] = ±2S ± 
\_S-,S+] = S z 



commutation relations for SL(2) triple (64) 



where SL(2) is a Lie algebra (4) . 

We wish to define the total angular momentum operator S 2 in terms of the angular momenta operators around each axis. It 
will be convenient to define the vector of angular momenta: S = (S x , S y , S z ) where S x = S++ S- and S y = i(S + - We note 
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that S 2 = S-S = T,ae{x,y,z} satisfies [S 2 , S z ] = 0. Since S 2 and S z commute and are both Hermitian, they are simultaneously 
diagonalizable, i.e., they share a common orthonormal eigenbasis. 

Recalling some basic results from the quantum theory of angular momenta, we note that for the basis {\S, mj)} where S 
represents the total spin quantum number of a set of spin-1/2 particles and ms represents the total spin projection quantum 
number onto the z-axis, we can show that 

S 2 \S,ms) = S(S+l)\S,m s ) (65a) 
S z \S,m s ) = m s \S,m s ) (65b) 

where 

56 H'4 2 '-"'y} (66) 

and 

m s e{-S,-S+l,...,S-l,S}. (67) 

Keeping in mind that the basis states of the good subspace are eigenvectors of the interaction Hamiltonian and also satisfy 
Eq. ( |32| ) for a e {+, -, z}, let us examine the cases N = 1,2,3 and 4 in turn. 



A. One Physical Qubit 

For a single physical qubit (N = 1), the basis {|0) , |1)} corresponds to that of our familiar spin-^ particle, with S = | and 
ms = ±|. We identify our logical zero and one states as follows 

|0) = |5=i m s =\) (68a) 
\1) = \S=\,ms = -\) (68b) 



B. Two Physical Qubits 

For two physical qubits (N = 2), which we label A and B, with individual spins Sa = h an d Sb = |, we first note that the 

prescription for adding angular momentum (or spin) given Sa and Sb, is to form the new spin operator S = Sa + Sb with 
eigenvalues 

S6{\Sa-Sb\,...,Sa + S b } (69) 
with the corresponding spin projection eigenvalues 

m s e{-S,...,S}. (70) 

Thus, for two physical qubits, we see that the total spin eigenvalues 5(jv=2) can only take the value or 1. For 5(jv=2) = 0, 
we see that mg can only take the value (singlet subspace) whereas when S^ =2 ) = 1, ms can take any one of the three values 
-1,0,1 (triplet subspace). For our singlet subspace, 

\S (N=2) = 0, m s = 0) = (|01) - |10)) (71) 

we see that S z \S^ =2 ) = 0; m S = 0) = for our system operator S z = a\ ® I + I ® a\. In fact, S a \S^ =2 ) = 0, ms = 0) = for 
a e {+, -, z) where S a = Y,i Similarly, it can also be shown that S 2 |5 , (jv=2) = 0, ms = 0) = 0. Since the singlet state clearly 
satisfies condition ( |32) , we conclude that |S'(jv=2) = 0, ms = 0) is by itself a one-dimensional DFS. However, we also note that 
the triplet states are not eigenstates of S Z ,S + , S- and thus violate Eq. ( |32| ). 



12 



C. Three Physical Qubits 



For three physical qubits (N = 3), let us label the physical qubits A, B and C each with corresponding total spins Sa = \, 
Sb = \ and Sc = \ ■ If we think of this system as a combination of a pair of spins (A and B) with another spin C, we can again 
apply our rule for adding angular momenta which gives us from combining our pair of physical qubits into a S(jv=2) = system 
with a spin- 1 particle, eigenvalues of the total spin operator of 



S, 



(N=3) 









1 


o-i 




o + i 




2 


, . . . , 


2 


~ 2 



with corresponding spin projection eigenvalues ms = ±|. If instead we chose to combine our pair of physical qubits A and B 



into a S, 



(N=2) 



1 system with a spin-^ particle, the eigenvalues of the total spin operator would be 



S, 



(JV=3) 



with corresponding spin projection eigenvalues ms 



1 3 

2' 2 

for S(jv=3) = \orm s = ±|,±i 



1 




1 


1 - - 




1 + - 


2 


J * ' ' ) 


2 



forS* ( 



(JV=3) 



These distinct cases 



arise because there are 2 distinct ways we can get a total spin of S = | from a system with 3 physical qubits, either with two of 
the qubits combined as a spin-1 system and then combined with the spin-| particle or alternatively with two qubits combined as 
a spin-0 system and subsequently combined with the remaining spin-| particle. 



D. Generalization to N physical qubits 

The extension of this idea of combining spin angular momenta is straightforward. There is an inductive method of building 
up from the above procedure to higher N. Suppose we wish to build up the spin states of TV physical qubits. We would first 
build up the states for a set of TV - 1 physical qubits and then couple the spin of the last qubit. 

Suppose we consider the case with N = 4 physical qubits. We can create a Bratelli diagram for this scenario (Figure |2j. The 
decoherence free states lie on the axis where S = 0. There are two possible paths to build up the states from N = to N = 4. So 
we can construct a qubit with each logical state |0l) and equal to a decoherence free state indexed by the path label A. 



\Q L )=/^\= \S=0,m s = 0,A = 0) 
\1l) = /\ = \S = 0,m s = 0, A = 1) 



(72) 



(73) 



Before proceeding, let us establish some notation. Define the singlet state |s)^ over the i and j qubit as 

1 



and the three triplet states as 



\t-) l ^\hl J ) = \S=l,m s = -l) 
1 



1*°), 



n/2 



(|0 i l i ) + |li0 i )) = |S = l,m s = 0) 



|* + > ii = l0 i J > = |5 = l,m s = l) 
For N = 4 physical qubits, the logical zero |0l) is given by 

= |singlet) ® |singlet) 

= l S )l2® l S >34' 

On the other hand the logical one will be later seen to be given by 



4= [IO12 ® l*~>34 + IOl2 ® l* + ) 3 4 " I*°)l2 ® l*°>3 4 ] 



(74) 

(75a) 
(75b) 
(75c) 



(76a) 
(76b) 
(76c) 



(77) 
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FIG. 2. Example Bratelli diagram for TV = 4 physical qubits. The decoherence free states lie on the points of the axis where S = 0. There are 
two ways of getting to S = when TV = 4 because there are two possible paths starting from TV = 0. So we can realize a qubit by setting the 
logical computational basis states |0l) and |1l) to these two decoherence free states. Parameter A indexes the two possible paths for TV = 4 
physical qubits. 



In a similar fashion, for TV = 6 physical qubits, we have for the logical zero 

|0i> = k> 12 ®|s} 34 ®|s> 56 (78) 

Note that permutations of qubit labels are permissible and can be used to define alternative basis states. Actually, we shall see 
in the next section that such permutations can be used to implement logical operations on the logical qubits. 



E. Higher Dimensions and Encoding Rate 



Clearly, more paths exist as TV grows. We thus have potentially more logical states available as we increase TV because the 
logical states correspond to the distinct paths leading to each decoherence free state. There exists a combinatorial formula for 
the number of paths to each point in the Bratelli diagram with S = for a given TV and hence for the dimension of the DFS 
T-L(N) of TV spin-| physical qubits 

TV! 

d N = dim (DFS ( TV) ) = - — — — ; rr (79) 

v v " (N/2)\(N/2 + l)\ 

As in the case of collective dephasing [Eq. ( [54] )] we can determine the encoding rate from the above formula. The encoding rate 
is the number of logical qubits Nl we obtain divided by the the number of physical qubits TV we put into the system. We can 
construct logical qubits from the logical states in the decoherence free subspace H(N), and the number TV^ of logical qubits is 
logarithmic in the number of logical states of T~L(N) with TV^ = log 2 (dAr). So the encoding rate rjy is 

_# of DFS qubits in H (TV) TV L log^ 

r/v = = = ■ (80) 

# of physical qubits TV TV 

It can be shown using Stirling's approximation 

log 2 TV! * (TV + 1/2) log 2 TV - N (81) 

for TV » 1, that the rate 

rN * l -2^r (82) 
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for N » 1 and hence that the rate rjy asymptotically approaches unity 

lim r N = 1. (83) 

This implies that when TV is very large, remarkably we get about as many logical qubits out of our system as physical qubits we 
put into the system. 



F. Logical Operations on the DFS of Four Qubits 

Consider once more the collective decoherence model. As we saw in Eq. (|62]> the systen>bath Hamiltonian has the form: 



H S b = £ S a ® B a 



where 



a=+,-,z 



N 



(84) 



(85) 



Recall that for N = 4, we have two DF states with 5 = 0. When we work out the states in the computational basis we find 
from Eqs. and $77\ that the corresponding logical qubits for this DFS are: 

|0) = - (|0101) - |0110) - |1001) + |1010» (86a) 



II 



) = (|1100) + |0011) - ||0101> - |lono> - 1|1001> - ||ioio>) 



(86b) 



Suppose that we have a string of information qubits that we wish to protect. One way we can achieve this goal is to implement 
a DFS (an error avoiding code). For our case, this is to group the qubits into blocks of length of 4, and encode each block into 
the logical qubits stated above. 

Now, Vx, y e {0, 1}, define the exchange operation Eij on the state \x)i ® \y)j to be: 



~ \v)i9\x)j 

Thus, has the following matrix representation in the standard basis of two qubits: 



(87) 



Ej 



I 1 \ 

10 

10 

V o o o i / 



(88) 



and it is trivial that 



[HsB,Eij]=Q 



(89) 



In the following we are going to see that E^j's can be used to generate the logical X,Y,Z operators for the logical qubits in 
our case. 

Consider the operator (-E 12 ), and observe its operation on the logical qubits: 



(~£i 2 )|0) 



- (E 12 ) X - (|0101) - |0110) - |1001) + |1010)) 
-i (|1001) - |1010) - |0101) + |0110)) 



10) 



- (E 12 ) -L (|1100> + |0011) - ||0101> - ||0110> - ||1001> - ||ioio>) 
~ (|1100) + |0011) - ||1001> - 1|1010> - ||0101> - ||0110>) 



-II) 



(90a) 

(90b) 
(90c) 
(90d) 

(90e) 
(90f) 
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We can see that (-£'12 ) acts like a Z operator on the logical qubits. 

Similarly, we may check that 4| (-E23 _ £-13) acts like an X operator on the logical qubits. 
Thus, we may define one set of X,Y,Z operators for the DFS in our case to be: 

a z = (-E 12 ) 

& x = (E23 ~ Ei 3 ) 

-V— ° [ -X -21 

a = 2^ 'J 



(91a) 
(91b) 

(91c) 



As we saw in Section III C with the <r x and a z operations we can construct arbitrary qubit rotations via the Euler angle 
formula: 



exp(i9h ■ a ) = exp(iaa x ) exp(iftcr z ) exp(i"fa x ) 



(92) 



To perform universal quantum computation we also need to construct entangling logical operations between the logical qubits. 
This too can be done using entirely using exchange operations. See [5 1 for the original construction of such a gate between the 
logical qubits of the 4-qubit DFS code. 

V. NOISELESS/DECOHERENCE FREE SUBSYSTEMS 
A. Representation theory of matrix algebras 

We begin this section by stating a theorem in representation theory of matrix algebras. 

Recall the general form of the system-bath Hamiltonian, Hsb = T, a S a ® B a . Let A = {poly(S' Q )} be the associative 
algebra generated by all the system operators S a (all sums and products of such operators). The system Hilbert space is 

H S = (C 2 )® N = C 2 ". 



Theorem 2. Assume that A is f -closed (i.e. a e A => e A) and that I e A. Then 

A=®I njXn] ®M dj xd 3 



(93) 



Here M.dxdi'C) denotes the algebra of complex-valued d x d matrices { M^xdiC) }, while as usual / is the identity matrix. The 
number j is the label of an irreducible representation (irrep) of A, rij is the degeneracy of the jth irrep, and dj is the dimension 
of the jth irrep. 

From Theorem [2] we know that every S a has a matrix representation in the form S a = ©j I nj x nj ® Md jX dj ■ 



S a 



Jl 











J2 











J 3 



(94) 



where J 7 - = I nj xn 3 ® M dj xd 3 - 



M d . xd . 
M djKdj 
M dj x d] 



(95) 



Therefore, we have the following Corollary resulting from the previous Theorem: 
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Corollary 1 (0). The system Hilbert space can be decomposed as 

^s = 0C™ 3 <s>C^, (96) 
i 

and tlie subsystem factors C nj 's are unaffected by decoherence. 

The above result can be easily seen since the system operators S a = ©j I nj xnj ® Mdjxdj act as the identity operator (hence 
do nothing) on each of the terms C" J . This result is important since we now have decoherence-free subsystems, which we can 
use to store protected quantum information. 

Recall that the "error algebra" A = {poly(S a )} acts trivially on the components C" J in Eq. |96|. Thus we can noiselessly 
store quantum information in one of the components C™ 7 . These components C" J are called noiseless subsystems (NS). Note 
that the dimension of Hs = (C 2 )® N = C 2 can be decomposed as 

2 N = Y i njd J . (97) 

J 

Define the commutant A' of A to be the set 

A' = {X :[X,A] = 0,VAeA}. (98) 

This set also forms a f-closed algebra and is reducible to, over the same basis as A, 

A' *®M njXnj (C)®I djxdj . (99) 
./ 

These are the logical operations for performing quantum computation: they act non-trivially on the noiseless subsystems C nj . 

B. Example: collective decoherence revisited 

1. General structure 

Let's return to the collective decoherence model. Recall that collective decoherence on N qubits is characterized by the system 
operators S a = £ i=1 erf , for a e { + , -, z}. In this case, the system space is 

N/2 

^s= C" J ®C dj , (100) 

J=0(l/2) 

where J labels the total spin, and the sum is from J = 0orJ = 1/2 if N is even or odd, respectively. For a fixed J, there are 
2 J + 1 different eigenvalues of raj, and hence 

dj = 2J+l. (101) 
By using angular momentum addition rules, one can prove that 

(2J+1)AH 

n ' J ~ (N/2 + 1 + J)\(N/2 - J)l ' ( 0Z) 

which is equal to the number of paths from the origin to the vertex (N, J) on the Bratteli diagram. For example, when N = 3: 

n i/2 = Mr = 2 - We have 

•H s = C™°<g>C 1 eC ni ®C 3 e- (103) 

for N even, and 

"H s = C ni/2 ®C 2 eC" 3/2 <g>C 4 e-", (104) 

for N odd. Note that when J = (so that dj = 1), C 1 is just a scalar and C n ° <8> C 1 = C"°, in which a DFS arises. 

The noiseless subsystems corresponding to different values of J for a given N can be computed by using the addition of 
angular momentum, as illustrated below. 
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2. The three qubit code for collective decoherence 
The smallest N which encodes one qubit in a noiseless subsystem is N = 3. In this case, 



H (N-3) = c 2 ®C 2 eC 1 ®C 4 , (105) 



Thus we can encode one qubit in the first factor C 2 of J = 1/2. The two paths of |0) and |1) are respectively (A = 0) 



and (A = 1). The end points of these two paths each have two spin projections mj = ±1/2 (since they correspond to 

a total spin J = 1/2). Using the state notation | J, A, raj), we thus have 

|0) =a\l/2, 0, -1/2) + [1/2, 0, 1/2) = |l/2, 0) ® (a |-l/2) + f3 |l/2», (106a) 
|1) =« |l/2, 1, -1/2) + (3 [1/2, 1, 1/2) = [1/2, 1) « (a |-l/2) + /3|l/2)), (106b) 

where a and /? are completely arbitrary. Or using the vector form, we have 

Suppose we want to encode a state \i}>) = a |0) + b |1). The encoded state is 

0) = a|6) + 0|l) = (j \ 9 l a V (108) 

where we only care about the encoded information a and b. Notice how this last result precisely corresponds to the C 2 ® C 2 term 
in Eq. ( |105| l. Thus, a and j3 are a "gauge freedom"; their precise values don't matter. 

The interaction Hamiltonian restricted to the system S is of the form 

3/2 

Hsb\ s = I nj ®M dj (109a) 

J=l/2 

=h ® A^ 2 ©It ® Ma (109b) 

(109c) 



I 2 ® M 2 






Mi 



What this means is that the term I 2 ® M2 acts on [i^) and leaves its first factor alone (this is good since that's where we store the 
qubit), but applies some arbitrary matrix M2 to the second factor (we don't care). M4 acts on the C 1 ® C 4 subspace, where we 
don't store any quantum information. 

We can check that the dimensions satisfy Eq. d97|: 



3/2 

E njdj = n 1/2 d 1/2 + n 3/2 d m = 2-2 + 1-4 = 8 = 2 3 . (110) 

J=l/2 

Let's now find explicit expressions for the basis state of the three-qubit noiseless subsystem. Recall that |0) = | J = 1/2, mj = 
1/2), |1) = [1/2,-1/2), the singlet state \s) = |0,0) = ^(|01) - |10)) and the triplet states are |t + ) = |1,1) = |00), |t_) = |1,1) = 

1 1 1 ) , and |to ) = |1j 0) = |(|01) + [10) ) . We now derive the four J = 1/2 states by using the addition of angular momentum and 
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Clebsch-Gordan coefficients. 



|l/2,0,~l/2)=| S )®|m 3 = -l/2) 
=-L(|011)-|101)) 
|1/2,0,1/2)=|»)®|0) 



V2 



(|010) -|100)) 



|l/2, 1, -1/2) =-^(V2| Jia = l,m Jl2 = -1) ® |m 3 = 1/2) - |J 12 = l,m Jl2 = 0) ® |m 3 = -1/2)) 



(2|110>-|011)-|101» 



|l/2, 1, 1/2) =^(ki2 = l,mj u = 0) ® \ma = 1/2) - V2\J 12 = l,m, 7l2 = 1) ® \m 3 = -1/2)) 



(|010) + |100> - 2 1001>). 



These are the basis states that appear in Eq. ( |106| >, so they complete the specification of the three-qubit code. 



(111a) 
(111b) 
(111c) 
(Hid) 

(llle) 

(lllf) 

(Hlg) 
(11 lh) 



3. Computation over the three-qubit code 

Consider the permutation operator E.^ = |(7 + tr, • dj) such that E^ \x) i |y) ■ = for x, y e {0, 1}. We have 

£ 12 |l/2,0, -1/2) =-L(- |011) + |101» = - |l/2,0, -1/2) (112a) 

E 12 |l/2, 0, 1/2) =4=(|100) - |010)) = - \l/2, 0, 1/2) (112b) 

E 12 |l/2, 1,-1/2) =^(21110)- |011) -|101)) = |l/2,l,-l/2) (112c) 

£fia |l/2, 1, 1/2) =^=(|010) + |100) - 2 1001)) = |l/2, 1, 1/2) . (112d) 
V6 

Thus E12 works as a logical -a 2 ', in the sense that 



E 
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1 

1_ 

in the ordered basis of the four J = 1/2 states given in Eq. ( |112| >. Again, this agrees with the C 2 
subspace where we store our qubit. 
Similarly, one can easily verify that 



(113) 



structure of the Hilbert 



1 

71 



Then a v can be obtained from 



(-B13--E23) = o- x ® I = a x 



2ia y = [d z ,a x ]. 



(114) 



(115) 



C. Intermediate Summary 



To conclude, the idea of a DFS/NS is powerful: we can use naturally available symmetries to encode and hide quantum 
information, and we can compute over the encoded, hidden information. But often such symmetries are imperfect, and we 
need additional tools to protect quantum information. Such an approach, which adds active intervention to the passive DFS/NS 
approach, is dynamical decoupling. 
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VI. DYNAMICAL DECOUPLING 



As we saw in the discussion of noiseless subsystems, the error algebra A = {poly(S' a )} is isomorphic to a direct sum of 
copies of dj x dj complex matrix algebras M.d,^d, (C): 

A = {Bl njXnj ®M djXdj (C), (116) 
j 

where nj is the degeneracy of the Jth irrep and dj is the dimension of the Jth irrep. We can store quantum information in a 
factor C" •' when nj > 1. However, from general principles (Noether's theorem) we know that degeneracy requires a symmetry, 
and in our case we would only have nj > 1 when the system-bath coupling has some symmetry. When there's no symmetry at 
all, nj = 1 for all J's, and a DFS or NS may not exist. Starting in this section, we discuss how to "engineer" the system-bath 
coupling to have some symmetry. 



A. Decoupling single qubit pure dephasing 

1. The ideal pulse case 

Consider a single qubit system with the pure dephasing system-bath coupling Hamiltonian 

H SB = <T Z ®B Z (117) 

and system Hamiltonian 

H s = X(t)a x . (118) 

We assume that X(t) is a fully controllable field, e.g., several pulses of a magnetic or electric field applied to the system. Assume 
these ideal pulses last for a period of time S, and with strength A, and 

SX=^. (119) 

Assume that at t = 0, we turn on the pulse for a period of time S, then let the system and bath interact for a period of time r, and 
repeat this procedure, as shown in Fig. [3] In the ideal case, 5^0 and A -*■ oo while still satisfying <5A = ^, which means the 
pulses are a series of delta functions. For simplicity we temporarily assume that H B = 0. 

X ♦ 



FIG. 3. Schematic of a dynamical decoupling pulse sequence. Pulses have width S and intervals of duration r. The modulation function \(t) 
is responsible for switching the pulses on and off. 

To formalize this "ideal pulse" scenario, let us define the system-bath "pulse-free" evolution operator f T , and the unitary 
transformation caused by the pulse, X, as follows: 

f T = e- iTHsB . (120) 

X = e'^"" ®I B = -ie~ lf,Xa * ®I B = ~ic x ® I B - (121) 

In the case of an ideal pulse (5 -*■ 0, A -*■ oo) there is no system-bath interaction during the time the pulse is turned on, since 
the duration of the pulse is 0. Then the joint system-bath evolution operator at time t = 2r is (dropping overall factors of i and 
minus signs) 

Xf T Xf T = a x e-" HsB a x e- iTHsB = e^™' ' HsB °* e -' lTHsB . (122) 
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where in the second equality we used the identity 



Ue A U f = e UAU \ (123) 



valid for any operator A and unitary U. 

On the other hand, since the Pauli matrices are Hermitian and every pair of distinct Pauli matrices anticommutes, 

{<J a ,vp} = 0, a* (3, (124) 

where the anti -commutator is defined as 

{A, B} = AB + BA (125) 
for any pair of operators A and B, it follows that the sign of Hsb is flipped: 

(J x H SB a x = a x a z a x ® B z (126a) 

= -a z ®B z , (126b) 

= -Hsb- (126c) 

This means that the evolution under Hsb nas been effectively time-reversed! 
Indeed, if we now substitute Eq. ( |126c| i into Eq. \\22) we obtain 

Xf T Xf T = e +iTH SB e -iTH SB = J (127 ) 

Thus, the bath has no effect on the system at the instant t = 2r. In other words, for a fleeting instant, at t = 2r, the the system is 
completely decoupled from the bath. Clearly, if we were to repeat Eq. ( |122| i over and over, the system would "stroboscopically" 
decouple from the bath every 2r. 



2. The real pulse case 

Unfortunately, in the real world, pulses cannot be described by 5 functions, because that would require infinite energy. Gen- 
erally, the pulse must be described by some continuous function X(t) in the time domain, which may or may not be a pulse (see, 
e.g., Ref. Q, for a variational solution of the optimal modulation for the pure dephasing problem). Then, during the period when 
the pulse is applied to the system, the system-bath Hamiltonian cannot be neglected, so we must take it into account. Keeping 
the assumption Hb = for the time being, we have to modify the pulse to 

X = g-^Ao-'+Hss) (128 ) 

If A » \\Hsb || and 5\ = ir/2 (we'll define the norm momentarily), then it's true that X « a x <S> Is, i.e., we can approximate the 
ideal pulse case of Eq. ( |121| >. Let us now see how good of an approximation this is. 

To deal with the real pulse case, we first recall the Baker-Campell-Hausdorff (BCH) formula (see any advanced book on 
matrices, e.g., IS]): 

e <A+B) = e eA e eB e (e 2 l2)[A,B]+0(e 3 ) ^ (129) 

for any pair of operators A and B. Now, set e = -id, A = \<r x , B = Hsb = cr z B z . Then the real pulse is 

X = e -iS^ x e -iSH SB e -8 2 \[cr*,HsB]/2+0(6 3 ) (13Q) 
ideal pulse OK does damage 

The first exponential term is just the ideal pulse, and the second is OK as well (we will see that shortly), but the third term will 
cause the pulse sequence to operate imperfectly. Let's analyze the pulse sequence subject to this structure of the real pulse. 
First, let us define the operator norm: 

H A H ESup TTiT\ ' (131) 

IV) vPw) 

which reduces to the absolute value of the largest eigenvalue of A when A is Hermitian. The operator norm is an example of a 
unitarily invariant (ui) norm: If U and V are unitary, and A is some operator, a norm is said to be unitarily invariant if 

\\UAV\\ m =\\A\\ m . (132) 
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Such norms are submultiplicative over products and distributive over tensor products [9|: 

||AB|| UI < ll^l|ui||B|| UI , \\A®B\\ ui =\\A\\ ui \\B\\ ui . (133) 
Then, using \\cr x \\ = 1 (the eigenvalues of a x are ±1), we have 

\\H SB \\ = \\o x ®B Z \\ = \\B Z \\. (134) 

Using this we find 

\\S 2 X[<j x ,H SB ]/2\\<6 2 X(\\a x H SB \\ + \\H SB a x \\)/2 (135a) 
<(tt/2)8\\(t x \\\\Hsb\\ (135b) 
= 0(<5||5 2 ||), (135c) 

where we used the triangle inequality and 5X = ir/2. So, we arrive at the important conclusion that the pulse width should 
be small compared to the inverse of the system-bath coupling strength, i.e., 6 « l/||i? z ||, should be satisfied assuming ||-B Z || 
is finite. This assumption won't always be satisfied (e.g., it does not hold for the spin-boson model), in which case different 
analysis techniques are required. In particular, operators norms will have to be replaced by correlation functions, which remain 
finite even when operator norms are formally infinite (see, e.g., ifTOl ). But, for now we shall simply assume that all operators 
norms we shall encounter are indeed finite. 

Let's Taylor expand the "damage" term to lowest order: 

e -s 2 \[**,HsB]/2 = I _ s 2 \[a x ,H SB ]/2 + 0(S 3 ) (136a) 
= I + 0(5\\B Z \\). (136b) 

Putting everything together, including e ~ l5Xa = -icr x , the evolution subject to the real pulse is, from Eq. ( |130| > (again dropping 
overall phase factors) 

Xf T Xf T = [o-^e-^^^CJ + 0(5|| J B 2; ||))]e- i7 " jHrss [cr 2: e-^ jH " ss (/ + 0(^11 1| )) ]e"^ jH " ss (137a) 

= e -i(T + &)a*HsB<y* e -r(T+8)H SB + 0(S\\B Z \\) + 0(S 2 \\B Z f) (137b) 
= ^(r^Hss^ir^HsB + 0(8\\B Z \\) (137c) 
= I + 0(S\\B z \\), (137d) 

so we see that the real pulse sequence has a first order pulse width correction. 

Now let us recall that in fact H B + 0. How does this impact the analysis? Both the free evolution and the pulse actually 
include H B : 

f T = e -"(H S B + H B )^ (138a) 

so we need A » || H$b + H B \\. Set H' SB = H B + H$ B , and note that the ideal pulse commutes with H B , so that 

a x (H SB + H B )a x = -H SB + H B . (139) 
Substituting Eqs. ( |138a[ > and ( |138b| i into Eq. ( |137[ ) we then have: 

Xf T Xf T = e ~i(-r+5)<7*H' SB ** ■ e -i{r+5)H' SB + o(6\\H' SB \\) (140a) 

= e -i(T+6)(.-H SB +H B ) e -i(T + 6)(HsB+H B ) + Q^jJJ^ ( 140b ) 

Setting A = H SB + H B , B = -H SB + H B , and using the BCH formula j!29\ again, we have A + B = 2H B and || [A, B] \\/2 < 
\\H B ~ H SB \\ \\H b + H SB \\ < (\\H S b\\ + l|gi?ll) 2 , so that Eq. ggg reduces to 

Xf T Xf T = I s ®e- 2l ^ HB +0[{t + 5) 2 {\\Hsb\\ + \\H b \\) 2 ] + 0[5{\\Hsb\\ + \\H b \\)1 (141) 

Assuming that the pulses are very narrow, i.e., 5 « t (recall that we anyhow need this for ideal pulses), we can neglect 6 relative 
to t in the second term, and so the smallness conditions are 

S«t«1/(\\B z \\ + \\H b \\), (142) 

which replaces the earlier 6 « 1/\\B Z \\ condition we derived when we ignored H B . 
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B. Decoupling single qubit general decoherence 

Let us now consider the most general 1-qubit system-bath coupling Hamiltonian 

H SB = £ a a ®B a . 



Using the anticommutation condition Eq. ( |124| i we have 

o x H SB a x = a x ®B X -a y ® B y - a z ® B z 



(143) 



(144) 



so that the Xf T Xf T pulse sequence should cancel both the y and z contributions. The remaining problem is how to deal with 
the a x term in Hsb- 



Pi 



z 



X 



X 



t 



FIG. 4. Schematic of the pulse sequence used to suppress general single-qubit decoherence. This pulse sequence is sometimes called XY-4, or 
the universal decoupling sequence. 

Let us assume that the pulses are ideal (8 = 0). We can remove the remaining a x term by inserting the sequence for pure 
dephasing into a second pulse sequence, designed to remove the a x term. This kind of recursive construction is very powerful, 



and we will see it again in Section IX 
Let the free evolution again be 

f T = e - iTHsB . 

Then, after applying an X-type sequence, 

Xf 2T = f T Xf T = e -Vr(<r*®B* + H B ) + 

To remove the remaining a x ® B x we can apply a y-type sequence to f 2r \ 

= YXf T Xf T YXf T Xf T 
= Z f T X f T Z f T X f T . 

where as usual we dropped overall phase factors. Clearly, 



(145) 

(146) 

(147a) 
(147b) 
(147c) 

(148) 



so that at t = 4r the system is completely decoupled from the bath. This pulse sequence is shown in Fig.|4j and is the "universal 
decoupling sequence" (for a single qubit), since it removes a general system-bath interaction. 



VII. DYNAMICAL DECOUPLING AS SYMMETRIZATION 



We saw in Eq. ( |147| i that the sequence Z f T X f T Z f T X f T decouples a single qubit from an arbitrary bath (to first order). We 
constructed this sequence using a recursive scheme. In this section we would like to adopt a different perspective, which will 
help us generalize the theory beyond the single qubit case. This perspective is based on symmetrization ifTTl . 
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Notice that, up to a global phase, 



Zf T Xf T Zf T Xf T = (Zf T Z) (Yf T Y) (Xf r X) (If T I) . 



(149) 



On the right hand side of ( |149| l we see a clear structure: we are "cycling" over the group formed by the elements {/, X, Y, Z}. 
Note that because we are not concerned with global phases, this is not the Pauli groupFl Rather, the four element group is the 
abelian Klein group, whose multiplication table is given by 



X 


I 


X 


Y 


z 


I 


I 


X 


Y 


z 


X 


X 


I 


Z 


Y 


Y 


Y 


z 


I 


X 


Z 


Z 


Y 


X 


I 



Returning to the decoupling discussion, to see why the sequence in Eq. ( |149| l works, note that if we let A a = a a ® B a we 
have 



-iT(A*+A v +A*+H B ) 



Ifrl = fr = e 
Xf T X = e 

Yf T Y = e - iTr7VHay = e -iT(-A*+AV-A z +H B ) 

Zf T Z = e 



-ira Ha x 



-ira Ha" 



-ir{A"-A y -A z +H B ) 



-ir(-A x -A v +A~ +H B ) 



(150) 
(151) 
(152) 
(153) 



Using the BCH expansion (|129| l again, we see that when we add all four of the exponents they cancel all A a terms perfectly, so 
that the right hand side of (|149|l is just 



(Zf T Z) (Yf T Y) (Xf T X) (If T I) = B 



-HtH b 



(154) 



just like in Eq. ( |148| ). This is the first order decoupling we were looking for. 

From the right hand side of equation ( |149| l, we also gain some intuition as to what our strategy should be beyond the single 
qubit case. Again, define 



fr = cxp[-iT(H SB +H B )], 



(155) 



where now Hsb and Hb are completely general system-bath and pure-bath operators. Generalizing from Eq. ( |149| l, consider a 
group 

G = {9o,---,9k} (156) 
(with g = I) of unitary transformations gj acting purely on the system, and apply the following symmetrization sequence: 



U{T) = f\g]f T g 3 

j=o 

K 



FT e - iT {9) H SB9] + H B ) 
3=0 

e -iT( J ^X« glH 3Bgi+ H B ) +0 ( T 2^ (157c) 

where we used Eq. ( |123| l in the second equality and the BCH formula in the third. The total time taken by the sequence is 

T=(K+1)t. (158) 

Now, define the effective, or average Hamiltonian 



(157a) 
(157b) 



1 



K 



H 'sB=-^— A Y J 9]H SB 9j. 

3=0 



(159) 



Our strategy for general first order decoupling could be one of following: 



1 The Pauli group is the 16-element group {±1, ±X, ±Y, ±Z, ±il, ±iX, ±iY, ±iZ}). 
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1. Pick a group Q such that H' SB = 0. 

2. Pick a group Q such that H' SB = Is ® B' . 

The first of these is precisely what we saw for decoupling a single qubit using the Pauli (or Klein) group, i.e., Eq. ( |154| i. To 
see when we can achieve the second strategy (which obviously included the first as a special case with B' = 0), note that H' SB 
belongs to the centralizer of the group Q, i.e. 

H' SB e Z(Q) = {A\Ag = gA V 5 e5}. 

To prove this we only need to show that g^ H' SB g = H' SB for all g e Q, since this immediately implies that [H' SB , g] = V g € £7. 
Indeed, 

g ] 'H SB g = — L- f g1g]Hg B g jg (160a) 

= ^TTE(5 J 5) 1 " J ff5Bfe.9) (160b) 
= H' SB , (160c) 



since by group closure {gjg}f =0 also covers all of Q. 

The fact that H' SB commutes with everything in Q means that we can apply Schur's Lemma [4|: 

Lemma 1 (Schur's Lemma). Let Q = {gi} be a group. Let T(G) be an irreducible d- dimensional representation of Q {i.e. not 
all of the T(<7j) are similar to a block-diagonal matrix). If there is a dx d matrix A such that [A, g{\ = V gi e Q, then A oc I. 

Thus, it follows from this lemma that, provided we pick Q so that its matrix representation over the relevant system Hilbert 
space is irreducible, then indeed H' SB oc I s . 

For example, Hs = (C 2 ) 8 ™ = C 2 for n qubits; the dimension of the irrep should then be 2" in this case. Which decoupling 

group has a 2"-dimensional irrep over (C 2 ) 8 ™? An example is the rt-fold tensor product of the Pauli group: Q = {/<S>---<S>/, X® 
I (8) ••• ® 7, . . . , Z <8> ••• ® Z}. And indeed, this decoupling group suffices to decouple the most general systen>bath Hamiltonian 
in the case of n qubits: 

H SB = Y l (7 l 1 ®"®< 7 n"® Ba > ( 161 > 

a 

where a = {ai, . . . , a n }, and a, e {0, x, y, z}, with the convention that a = I. Fortunately, such a system-bath interaction is 
completely unrealistic, since it involves n + 14oody interactions. 'Fortunately," since the decoupling group we just wrote down 
has K - 1 = 4™ elements, so that the time it would take to apply just once symmetrization sequence ( |157a| > grows exponentially 
with the number of qubits. 

Actually, this approach using Schur's lemma is a bit too blunt. As we shall see, it is possible to drastically reduce the required 
resources for decoupling, for example by combining decoupling with DFS encoding, or by focusing on more reasonable models 
of system-bath interactions. 

VIII. COMBINING DYNAMICAL DECOUPLING WITH DFS 

We saw that to decouple the general system-bath interaction H$b in Eq. ( |161| l would require a group with an exponentially 
large number of elements. This is not only impractical, it might also destroy any benefit we would hope to get from efficient 
quantum algorithms. Therefore we now consider ways to shorten the decoupling sequence. As we'll see, this is possible, at the 
expense of of using more qubits. There will thus be a space-time tradeoff. For an entry into the original literature on this topic 
see Ref. US. 

A. Dephasing on two qubits: a hybrid DFS-DD approach 



Consider a system consisting of two qubits that are coupled to a bath by the dephasing interaction 

Hsb = <J 7 i®B{ + a! ® 2?2- 



(162) 
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This Hamiltonian is not invariant under swapping the two qubits since they couple to different bath operators. To make this more 
apparent, rewrite the interaction as 

# SB = (^)®£_ + (^)®S + (163) 

where the redefined bath operators are B ± = B\ ± Bf. We find that ) is a "collective dephasing" operator that applies 

the same dephasing to both qubits, while the "differential dephasing" operator ( gl ^ 2 j applies opposite dephasing to the two 

qubits. We now notice that the states |01) and 1 10) (and all superpositions thereof) are annihilated by the collective dephasing 
operator: 

(ZI + IZ)(a\01) + /3|10)) = V a,/3eC, (164) 

where ZI = o\ ® I and IZ = I <g> erf- These can therefore be chosen as basis elements of a DFS in which a single logical qubit 
can be encoded as |0) = |01) and |1) = |10), just as inEq. |47]). 

Now that we have chosen a basis that vanishes under the effect of one part of the interaction hamiltonian, this effectively 
reduces the interaction to 

ZI — IZ 

H S b\dfs = ^ ® B - = ° Z ® B - (165) 

If the initial interaction had been symmetric, choosing the DFS would have reduced it to zero. However, the interaction was 
not symmetric in this case, and we are left with the above differential dephasing term. We notice further that the residual term is 
the same as a a z , or logical Z operating on the DFS basis|^] We recall that dephasing acting on a single qubit was decoupled by 
pulses that implemented the X orY operators, and hence expect that the a z interaction can be decoupled using alorF pulse. 
We'll use the convention that logical/encoded terms in the Hamiltonian are denoted by a a , while the corresponding unitaries are 
denoted by X, Y, or Z. Thus 

,* = XX + YY ,v, YX ~ XY . (166 ) 
2 2 

Restricted to the DFS, the implementation of an X pulse using a a x is analogous to the implementation of an X pulse by 
applying a x for an appropriate period of time. 

e -if9- =e -i^(XX + YY) (16?a) 

(using [XX, YY] = 0) (167b) 



\XX -%%YY 



= e 4 e 

= -^=[I-iXX] — [I-iYY] (using (XX) 2 = (YY) 2 = I) (167c) 
v 2 v 2 

= - [II - iXX - iYY + ZZ] (167d) 

= ^(XX + YY) = -iX, (167e) 

where the term // + ZZ in Eq. ( |167d[ ) was ignored since it vanishes on the DFS. 
Hence, the dynamical decoupling process is effective in the sense that 

Xf T Xf T \ DF5 = exp[-2zr(7® B)] + 0[(2r) 2 ], (168) 

where B is a bath operator whose exact form does not matter, since we have obtained a pure-bath operator up to a time 0(T 2 ), 
where T = 2r. What have we learned from this example? That we don't need to remove every term in the system-bath 
Hamiltonian; instead we can use a DFS encoding along with DD. Next we'll see how this can save us some pulse resources. 



B. General decoherence on two qubits: a hybrid DFS-DD approach 

We now consider the most general system-bath Hamiltonian on two qubits: 

Hsb= £ K 1 ®^ 2 )®^ 1 " 2 , (169) 



2 This is a symmetrized version of the logical Z operator in Eq. (57) . 
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Effect on DFS states 


Operators 


unchanged 


II, Z I + IZ,ZZ,XX - YY,XY + YX 


logical op. 


— z —y — x 

a ,a ,<j 


leakage 


XIJX,YIJY,XZ,ZX,YZ,ZY 



TABLE I. Classification of all two-qubit error operators on the DFS for collective dephasing. 



where a,; e {0, x, y, z}. 

Within the framework of the same DFS as earlier (DFS=span{|01), |10}}), we can classify all possible (4 2 = 16) system 
operators as either 

• leaving system states unchanged (i.e., acting as proportional to /) 

• mapping system states to other states within the DFS; these correspond to logical operations (these are errors since they 
occur as a result of interaction with the bath) 

• transitions from the DFS to outside and vice versa ("leakage") 
The operators causing these errors and their effects are given in Table [j] 



For example, ZZ acts on |0) = |01) and |1) = |10) as while XI takes both |0) and |1) out of the DFS, to |11) and |00), 
respectively. 

Along the same lines as single qubit dynamical decoupling, we look for an operator that anticommutes with all the leakage 
operators, and an operator that anticommutes with all the logical operators. Performing a calculation very similar to Eq. \\61\ 
we find that 

• exp(-ma x ) = -ZZ and anticommutes with the entire leakage set 

• exp(-z^(T 2 ) = -iZ and anticommutes with the logical error operators a x ,a y . 

A combination of these operators, along with X which we used above, is sufficient to reduce the effect of the system-bath 
interaction Hamiltonian to that of a pure bath operator that acts trivially on the system. First we apply X to decouple the logical 
error operators a z and a y , giving us a net unitary evolution 

U 1 (2T)=Xf T Xf T + 0[(2 T ) 2 ]. (170) 

We already know that U\ does not contain a z and a y to first order in r, leaving us to compensate for the leakage errors. This is 
accomplished by using a ZZ pulse: 

J7 2 (4r) = ZZ-Ux{2t) -ZZ-U^t) +0[(4t) 2 ]. (171) 

All that remains now is to remove the logical error operator a x , since it commutes with both the X and ZZ pulses we have used 
so far. This can be performed using Z, which anticommutes with a x . Hence, the overall time evolution that compensates for all 
possible (logical and leakage) errors is of period 8r and is of the form 

U 3 (8t) = ZU 2 (4t)ZU 2 (4t) = e - %lT(1 ® B) + 0[(8t) 2 ]. (172) 

A qubit encoded into the |0) = |01) and |T) = 1 10) DFS is acted on (at time T = 8r) only by the innocuous operators in the first 
row of Tabled As a result it is completely free of decoherence, up to errors appearing to 0(T 2 ), while we used a pulse sequence 
that has length 8r, shorter by a factor of 2 compared to the sequence we would have had to use without the DFS encoding (the 
full two-qubit Pauli group). This, then, illustrates the space-time tradeoff between using full DD without DFS encoding, vs using 
a hybrid approach, where we use up twice the number of qubits, but gain a factor of two in time. 
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IX. CONCATENATED DYNAMICAL DECOUPLING: REMOVING ERRORS OF HIGHER ORDER IN T 



The dynamical decoupling techniques considered so far have all involved elimination of decoherence up to first order in time. 
We now consider the question of whether it is possible to further refine these techniques and remove the effect of noise up 
to higher orders in time. We saw in the earlier sections that applying pulses corresponding to the chosen decoupling group 
effectively causes a net unitary evolution 

k 

Udd{T) = U(d U f9i) = 9lU f 9o9lU f gi....glU f g k (173) 

where Uf is the free unitary evolution opreator Uf = e~ lHT °, g\ Uf 9 i = e~ lg ^ H 9iT ° , and T = kr . Without loss of generality, we 
rewrite H as 

H = H C + H NC (174) 

where H c commutes with the group Q and H NC does not. This allows us to simplify the evaluation of Udd- We proceed by 
using the Baker-Campbell-Hausdorff formula, which yields 

U^ = UglUfg t = e^UT Y / glHg l+ ^Zl9lHg^g}Hg 1 ] + O(T^\ (175) 

i [ i Z i<j J 

The first term in the BCH series (0(t o ) term) can be evaluated as 

£ g\Hgi = £ g\n c g l + £ g\n Nc9i = ( k + l)H c + Y l iH Nc9l (176) 

i i i i 

We observe that the term in the above equation lies in the centralizer of Q: 

gt(T,9i H Ncg t )gs = Y,9Wi H Ncg s gi = £sj#jvc& (177) 

i i j 

We therefore set = (k + l)Hc + Z g^H^cgj, which has the property [Hq\q] = 0, i.e., the first order correction in the 
BCH series commutes with all the decoupling group elements. With this redefinition, we can rewrite 

U™(T) = exp{-*i#Vo - W^t* + ...} (178) 

where H^ c comprises the higher order terms (at least 0(t 2 )) of which, in general, do not commute with Q. We now 
begin the process of concatenation, which essentially amounts to removing the lowest order term in r of H^ c using a second 
iteration of the group action. This new sequence, U^ D will have an effective Hamiltonian, such that H^ c is proportional 
to C(tq ) or higher. This can be repeated, until the mth iteration where is proportional to C(r™ +1 ) or higher. 

(2) 

To see this consider U DD , which is constructed as follows: 

k 

U ( nl(kT) = U9iUi(T) gi (179a) 

k 

= ft g\ exp {-iff^Ti, - + ..) 9i (179b) 

= exp| -i^gj (H^ro + H { ^) 9l + £ [gl(H^r + H NC ^)g u gi{H^T Q + H NC ^) 9j ] + ... 1 



(179c) 



Looking for terms in the expansion above that do not commute with Q and are responsible for noise, we find that the lowest 
order term of this sort is 0(tq ). The part that does commute is labelled and we get 

U^ D (kT) = e -^ 2) ^+<^o+o(^)...} (18Q) 
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In a similar way, we can obtain U3 from U2 as, 



U { D %(k 2 T) = e A Hl C 3) ^<%r^O^)...}^ (181) 

and in general, the CDD process can be defined and described as 

U^(T {m) ) ■■= UglU^iT^- 1 ^ (182a) 

= e -i{^ m> ro+<™ ) r™ +1 + 0(r m+2 )...} j (lg2b) 

where the total duration of the sequence is given by 

r (m) = fc m T ^ (lg3) 

From the above formulation, it seems that we can eliminate the effects of noise on the time evolution of the system to any 
desired level of accuracy, merely by performing repeated iterations of the dynamical decoupling procedure. Note, however, 
that in principle the norm of may grow with m and thus is it not obvious at this point that increasing orders of CDD 

concatenation implies a better performance. Let us show that in fact there is an optimal level of concatenation^] 

Given the unitary evolution generated by J7^(T) = eT l ^ H c t o+h c o(t )) one can (j e fi ne an " error phase", 



cf>cDD(n)=T\\H^l (184) 
with T the total time of the sequence. This error phase can be bounded as lfl3l Eq. (46)] 

<Wn) = Tfltf^fl < T(2" 2 (/3r )"J), (185) 

where n is the degree of concatenation, 

J=\\H S b\\, 0=pTfll. ( 186 > 
and it was assumed that j3 « 1/T. It follows then that for an n-level CDD sequence we have the following two possibilities: 

• Assume a fixed T = k n TQ, i.e., tq can be made arbitrarily small. Then 0cdd(«) ^ T(2 n (/3T/k n ) n J), and thus the noise 
strength decreases monotonically as tq decreases as long as (/3ro(2/fc)™) < 1 Vn. This is in particular true if (3T(2/k) < 1. 
In practice, what this implies is that for a fixed T, provided f3T is small enough, the more concatenations of DD the better. 

• On the other hand, if r has a minimum physically achievable value, is there an optimal level of concatenation, i.e., a 
degree of CDD concatenation such that the error is minimized? What is the behavior of </>cdd("-) < fc™To(2 n (/3to)™ J) 
as n grows? We see that 

log^cDoW) < log(fc"r (2" 2 (/?T )" J)) 

= nlog(fc/?r ) + log(2)n 2 + log( Jr ). (187) 

Whenever kfirQ < 1, this expression has a minimum at the optimal concatenation level n opt = - l °2\o^) • This means that 
for n larger than n opt the CDD process loses its effectiveness, in contrast with the previous case. 

X. DYNAMICAL DECOUPLING AND REPRESENTATION THEORY 

In this section, our goal is to illustrate the connections between first-order DD and the result from representation theory, 
Theorem|2] which is a theorem of fundamental importance in the theory of quantum error correction. See Refs. ifTTl [Pfl [TBI for 
entries into the original literature on this topic. 



3 See Ref. (TO) for a rigorous analysis; here we adapt the more intuitive presentation in Ref. 1131 . 
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A. Information storage and computation under DD 



We define the group algebra CG of the group Q over the complex field as: CQ = { all linear combinations of the elements in 
Q, over C}. For example, if Q = {I, a x ,<r y , a z }, the Pauli group, then 



eg = 1 a a° c 

I a=0 7 x,y,z 



(188) 



Now consider the group algebra of our decoupling group, CG- Clearly it is a matrix algebra of dimension d x d where 
d = climbs ). We can always choose the group in such a way that CG is f-closed. Since every group includes the identity, we 
can thus invoke Theorem[2j using which we have: 



cg = @i nj ®M dj , 
j 



(189) 



where J is the irrep (irreducible representation) label, nj is the multiplicity of irrep J and dj is the dimension of the irrep 
labelled by J. The system Hilbert space Hs is correspondingly partitioned into a direct sum of product spaces, which we can 
write as 



J 



(190) 



Hence, every DD pulse that we apply to the system acts like identity on C™ 7 and as some non-trivial operation on C J . It's clear 
that if we store our quantum information in C nj , then the pulses do not affect it (of course the system-bath interaction can still 
affect information stored in C nj ). 

Now consider the commutant of the group algebra CG [see Eqs. ( |98| ) and (j99]»], 



CG' = {A\[A,CG] = 0} 

= (&M nj ®I dj . 

j 



(191a) 
(191b) 



We can immediately see from this definition and from Eq. ( |159[ ) that, by linearity, the effective system-bath Hamiltonian, H' SB , 
lies in the commutant of the group algebra, i.e., H' SB € C^'. It follows that H' SB can be represented as: 



H' SB = ®{H' SB ) nj ®I dj , 
j 



(192) 



i.e. the effective system-bath Hamiltonian has this block-diagonal representation, with blocks labeled by the irrep index J, each 
of dimension njdj, and where each non-trivial factor (H' SB ) nj is an nj x n,j matrix. It is clear that in this case, we can 
encode our quantum information in C dj , since H' SB will act as identity on it. Note that it doesn't necessarily matter that the 
pulses themselves have a non-trivial effect on Cdj '■ because we know everything about the pulses, they can compensated for by 
applying appropriate transformations. 

However, it is often desirable not to have to compensate for the action of the DD pulses. In that case we might we want to be 
able to store information in C" 7 rather than C dj . If so, clearly we need to somehow make the effect of H' SB trivial on C nj as 
well. Since H' SB is determined by our group of pulses, G, it boils down to choosing the appropriate set of pulses, i.e., picking 
the group G such that 



H' SB e CG' n CG 

= 0Aj/ n; <s> I dj , Xj e 



(193a) 
(193b) 



where Eq. ( |193b| i can be easily deduced by examining Eqs. ( |189| l and ( |191b| i. This means that if we pick an appropriate G, H' SB 
will have the form: 



' Xil ni d 1 



h 'sb 



A,/, 



(194) 



Suppose in the above: 
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• We have nj = 1, then there would only be a single J in the direct sum, which would imply that dj = d. This would mean 
that H' SB oc I d , which was the result we obtained by applying Schur's Lemma ([I]). Thus, that result was a special case 
resulting from this more general structure. 

• If Xj = for all J, then H' SB = and the system bath interaction is annihilated (to first order in T). This describes the 
kind of situation we obtained with the Klein group symmetrization [see Eq. ( [148) ]. 

So, now we have established that it is possible to protect information stored both in C nj and C dj . Information stored in the 
latter is protected since the effective system-bath interaction acts like identity on that space; though we'd have to compensate 
for the influence of the pulses we would apply. Storing in the former requires us to cleverly choose the pulse group Q such 
that the effective system-bath interaction becomes a block-diagonal matrix, with each block proportional to identity. Our choice 
of storage location might depend on which of maxj nj and max,; dj is greater, since that would provide us with a larger- 
dimensional space and hence more qubits. 

We now discuss some examples where we will use these techniques: 



B. Examples 

In all the following examples, our system is a set of 7V-qubits, i.e., 1-Lg = (C 2 )® N = C 2 . 



1. Example 1: Q = (SU(2))® N 

In this case the pulses are products of N arbitrary single-qubit unitaries, and the only operator which commutes with Q is 
the identity operator, i.e., H' SB e CQ' = C. This is the Schur's Lemma situation again. Thus, in this case we have the choice 
of storing information in either the left-hand (C nj ) or right-hand factor (C dj ), or both if we do not mind compensating for the 
action of the pulses on the right-hand factor. This case is interesting since it can be well approximated by picking the SU(2) 
rotation on each qubit at random |[T6l . 



2. Example 2: Q = Collective 5(7(2) 



In this case we have added to Example 1 the constraint that the same unitary matrix acts on every qubit. Then H' SB e CQ' = 
CSn, where Sn is the permutation group on TY-elements. Thus the effective system-bath Hamiltonian is not proportional to 
identity, and instead we have here the case where H' SB acts trivially only on the right-hand factors. Its action on the left-hand 
factors is to apply permutations. 

The collective-SU(2) group is generated by the sum of the Pauli matrices, i.e., o't}a=x.y,z- We already encountered these 
sums in the study of collective decoherence [see Eqs. |62} and (|63]l]. In this case, our pulse group has the same generating 
Hamiltonian as the system part of the system-bath interaction in collective decoherence: iJ™^ dec - = Y, a =x y z(T,i <J i') ® B a . 
Since the DD group acts as identity on the left-hand factors, and the DD group behaves like the collective decoherence operators, 
while the effective system-bath interaction H' SB acts like the exchange operators ( [87] > we encountered in our study of computa- 
tion over the DFS for collective decoherence, we see that the current situation is the reverse (or dual) of the situation back in the 
DFS case. In other words, we can invoke the machinery we developed then [see Eqs. ( |100| i-( fT02] i] but we should flip the role of 

dj and nj: H s = ©jio(i/2) C "' ® C<i '' where now 

nj = 2J+l, (195a) 

dj = (N/2 + 1 + J)l(N/2 - J)! ' (195b) 
We can always pick an irrep J so that dj > nj, and so we have here the same code rates as in the case of a DFS. 



3. Example 3: Q = S n 

In this case Q is the permutation group. As we saw in the last example, the permutation group is dual to collective SU(2). 
So in this case, H' SB e CQ' ="collective decoherence". In fact, since the permutation group can be obtained by swaps (or 
transpositions), we could as well take Q = {SWAPij}. And we know that, 

SWAP.j = e -<**i-**/4 (196) 
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which is generated by the Heisenberg interaction, so the decoupling group is implementable in physical systems (such as quantum 
dots) where the exchange interaction is controllable (see Ref. ifTTl for a discussion of how to use efficiently implement Q in this 
case). The Hilbert space again splits as %s = ©j2o(i/2) ® where now we have the irrep dimension and multiplicity 
formulas we encountered during the DFS study: 

(2J+1)N\ 

nj = -, — , ttt — , rr (197a) 

J (N/2 + 1 + J)\(N/2 - J)l 

dj = 2J+l. (197b) 

We have the option of encoding into the right-hand factor, where H' SB acts as identity, but then the space dimension is only 
2J+ 1. Alternatively, we can encode into the left-hand factors, where the effective system-bath interaction has non-trivial action, 
but it acts as collective decoherence, so that our DFS encoding will completely hide the quantum information from the action 
of H' SB . This has the significant advantage (over right-hand factor encoding) of providing us with a code space of dimension 
(A72+^ 2 /j) l !(jv/2-./)! • This example leads us to the interesting conclusion that in this case in fact Eq. ( |193b[ ) applies, i.e. the 
effective system-bath interaction acts trivially everywhere. 



4. Example 4: Linear System-Bath coupling 
We consider a system-bath interaction of the form, 

Hsb= E £<®^f- 

ot=x,y,z i 



(198) 



Each qubit in this case has its own bath. In this noise model, we don't consider bilinear terms in the system such as af a ^ ® B 
because this is a 3-body interaction which is typically much weaker in nature and also very hard to engineer. 
The decoupling group we select for this is, 



a/3 



g _ jj®JV j£®N y® n Z® N } 



(199) 



We choose N to be even. Therefore, Q becomes abelian. And from representation theory, we know that all the irreps of an 
abelian group are 1-dimensional (scalars, so dj = 1 V J) and the number of irreps is the order of the group. Here \Q\ = 4. The 
irreps are: 



J 








1 


1 1 


1 


1 


2 


1 1 


-1 


-1 


3 


1 -1 


1 


-1 


4 


1 -1 


-1 


1 



The group algebra is, 



CQ = @I nj ®M dj 
j 

4 

= C,]I 2 N--2, 
J=l 



(200a) 
(200b) 



where in Eq. ( |200b[ ) rij = 2 N ~ 2 because T,j n jd,j = 2" and all the n/s are equal (by use of a standard multiplicity formula 
from group theory: exercise for the reader), we have 4nj = 2 N , making nj = 2 N ~ 2 . 

Thus the group algebra can represented as a block diagonal matrix, each block being proportional to an identity of dimension 
2 N ~ 2 with proportionality constant cj, i.e., we have (N - 2) qubits in each block that will be unaffected by the system-bath 
interaction. 

Let us pick the first (trivial) irrep to encode our qubits into, i.e., {1,1.1,1}. In this irrep each pulse acts as 1, so we're looking 
for code-states which are "stabilized" by the group (each group element acts as identity). After a bit of thought it is clear that 
such states are of the form: 



i 

7^ 



(|r) + |f» 



(201) 
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where f = NOT(r) and r e {0, 1}^ is an even weight binary string of N bits, i.e., r contains an even number of l's. Then it is 
easy to see that the action of any member of the decoupling group Q, leaves \ip r ) unchanged, so indeed \ip r ) belongs to the trivial 
irrep. 

Why did we pick this decoupling group? Because it has some very interesting and useful features which we list: 
Feature 1: We can show that H' SB = 0. 



Proof. 



H' S B = ^T9}H'sB9j (202a) 

\y\j=o 



J [IH SB I + X® N H SB X® N + Y® N H SB Y® N + Z® N H SB Z® N ] 
^[ H SB + £ i a f ® B f ~ "i ® B J - *i ® B f 



(202b) 



of ® Bf + of ® Bj - of ® Bf 



- of ® Bf - oj (8) Bj + of <S> Bf ) ] (202c) 
= 0, (202d) 

where to arrive at Eq. ( |202c) > we used the properties of the Pauli group to do the multiplication. For example, 

X %N of X %N = (of ® . . . ® of ® . . . ® of)(of)(of ® . . . ® of ® . . . ® of) (203a) 
= -a-Y, (203b) 

and so on. □ 

Thus, this decoupling group eliminates the system-bath interaction completely, to first order. Recall that this means that we 
can use the left-hand factor to encode and store information. This allows us to perform computation using the commutant, which 
has non-trivial action on the left-hand factor. 

Feature 2: We can do computation on the decoherence-protected qubits. 

For N = 2, i.e., for 2 physical qubits, we have no logical qubits, as the only one state possible according to Eq. ( |201[ ) is 
|t/v) = ^7g(|00) + 1 1 1 ) ) - This agrees with the fact that, since each irrep provides us with N - 2 logical qubits, we have zero 
logical qubits for 2 physical qubits. 

Let us list the possible states in the case of N = 4, using Eq. (|201[> 



l^oooo) = ^(l° 00 °) + I 1111 )) = I 55 ) . ( 204a ) 

IVoou) = -^(l° 011 > + I 1100 )) = I 15 ) > < 204b ) 

hfoioi) = -^(10101) + |1010)) = |01) , (204c) 

l^ouo) = ^(l° n °) + I 1001 )) = 111) » ( 204d ) 

and since we have 4 orthonormal states, we can use them as 2 qubits. Again, this agrees with the fact that TV - 2 = 2 in this case. 
In the above equations, the logical states have a bar on top of their labels. 

How do we perform computations on these states? For that we use the commutant of the group, 

CQ' = @M nj ®I dj (205a) 

= © cjM 2 n-2 (205b) 
j=i 

where c,j are scalars. We can check that the commutant can be generated by {XiXj+i}^ 2 u {Zj+iZjy}^ 2 . For example, 
we don't require Y^Yj to be in the generating set, since YjYj = XiXjZiZj. Therefore, for N = 4, the generating set for the 
commutant becomes { X\ X^ , X\ X3 , Z2Z4, Z3Z4}. 
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Let's check the action of X\X 2 on our first logical state IV'oooo) = |00), 

X t X 2 |00) = (X 1 <8> X 2 <8> I <s> 7) |-^=(|0000) + |llll»j (206a) 

= — (|0011) + |1100)) = |10). (206b) 
v2 

We can similarly check that X1X2 |01) = |Tl). Therefore X1X2 is logical Pauli-X on the first logical qubit. We write this as 
X1X2 = Xi, where the bar denotes a logical operator. Similarly we can verify that: 

XiX^Xj, (207) 
Z J+1 Z N = Z r (208) 

And having obtained logical-X and logical-Z (note that they anti-commute, as they should), we can implement any one qubit 
gate we like using only two-body interactions. And moreover, these logical gates lie in the commutant of the decoupling group. 
Therefore, we can compute while at the same time applying DD. 

But, in order to perform universal quantum computation, we also need to be able to perform entangling operations on two 
qubits, for example the CNOT gate. Let's examine XiXj which is a logical entangling operation : 

XiXj = XxX M XxX j+l (209a) 
= X i+1 X j+1 , (209b) 

which is a physical 2-body interaction. So, we have managed to implement a logical entangling operation using only a two-body 
interaction. 

It is well known that if we can implement any Hamiltonian of the form, 

H S = Y,u?{t)*? +E"f(*K z + E E (210) 

i i ae{X,Z} i,j 

then we can perform universal quantum computation [ Q. Therefore, in our case, replacing the Pauli operators in Eq. ( |2 10] > with 
their logical counterparts and expanding them in terms of their decompositions ( |207| i in the physical qubit space, we obtain: 

N-2 N-2 

H S = E + E sf(*)*w + E E Mt)°M (2Ha) 

i=l i=l ae{X,Z} i,j 

N-2 N-2 

= E w «(*VKi + E^(«K t /N+ E E^+u+i(^?+i^+i (211b) 

i=l i=l ae{X,Z} i,j 

Remarkably, we have obtained decoupling and universal quantum computation on the logical qubits using only two-body inter- 
actions on the physical qubits. A similar approach has been studied numerically in the context of the 4-qubit DFS code, with 
gates protected by CDD, showing evidence for a highly robust set of universal gates [ 18 1. Related ideas apply in the context of 
adiabatic quantum computation |fT9l . 



XI. CONCLUSIONS 



This review has covered a selection of topics in the theory of decoherence-free subspaces, noiseless subsystems, and dynamical 
decoupling. We have seen how these tools allow one to hide information from the environment, and when this hiding is imperfect, 
how dynamical decoupling allows us to suppress the remaining residual decoherence. Moreover, we have shown explicitly how 
universal quantum computation is compatible with decoherence avoidance and suppression. 

Many important topics were left out in this brief review. For example, apart from CDD we did not address high-order 
decoupling methods, in particular schemes based on optimized pulse intervals I1201 - I231 . Nor did we address the filter function 
approach to DD ll24l . optimized continuous modulation Gl[25], or randomized decoupling, which is well suited to strongly 
time-dependent baths 11261 . It is important to stress that beyond decoherence avoidance and suppression, the theory of noiseless 
subsystems gave rise also to important advances in the theory of quantum error correcting codes, such as operator quantum 
error correction [27|. However, perhaps our greatest omission has been the abundance of experimental results which have both 
confirmed and driven the theoretical developments described here. Nevertheless, hopefully we have given the reader the tools 
and inspiration to delve deeper into the large and fascinating literature on decoherence avoidance and suppression. 
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